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ABSTRACT 

The process e + e~ — > u, + fi~ tfx, where the missing energy is due to the production 
of a tower of invisible graviton states in a model with large extra dimensions, is 
considered. We focus on the scenario when this process is used to detect a heavy 
dileptonic resonance in the electroweak sector of the model, taking as example 
some models with an extra Z' boson. It turns out that at a 500 GeV machine 
with 1000 ftr 1 of luminosity, it may be possible to use this process to detect such 
resonances if there are two large extra dimensions and the string scale is not too 
far above a TeV. 



A concerted international effort is currently under way to build a high energy linear e + e~ 
collider, which, in the wake of expected discoveries at the LHC, would conduct precision 
experiments to probe the TeV scale. While discussions on the technical design and feasibility 
of such a machine are already under way[l , it is essential, at the moment, to study the physics 
possibilities of such a facility and to determine what would be the best options to have when 
the experimental design is finalized. Quite naturally, the study of such physics possibilities [2 
is heavily dependent on the experience of the past, where the LEP collider was specifically 
designed to refine and complement the discoveries made with the proton synchrotron at CERN. 
The difference, in this case, is that we have no unique prediction of the discoveries expected at 
the LHC, and hence, the design for a linear collider must be broad-based and robust enough 
to adjust to all contingencies. 

Ever since the first serious accelerator experiments in the 1960s, resonance hunting has always 
been a major concern of any high energy experiment. A hadron collider is an ideal machine 
for seeking resonances, as the natural variations in parton momenta provide a wide spread 
of effective centre-of-mass energies, which then allow the machine to pick up any unknown 
resonances within its kinematic reach. This is not the case at an e + e~ collider, where the initial 
state energies are fixed and have to be decided as part of the machine design. This feature has 
obvious advantages for precision measurements, as, for example, at LEP-1, where the energy 
was fine-tuned to the Z° resonance, but there is always the possibility that there may exist an 
unknown resonant state which lies several decay widths away from the actual machine energy, 
and consequently, produces a very small effect. A major concern at a linear e + e~ collider, 
would, therefore, be the identification of ways and means of 'seeing' such resonances. This 
can be achieved if there is a mechanism to spread out the initial state energies. Fortunately, 
a mechanism is readily at hand, namely initial-state radiation and beamstrahlung, producing 
either hard, transverse photons (which can be tagged) or soft, collinear photons which will go 
unobserved down the beam pipe creating a mismatch in the observed momentum balance in 
the longitudinal direction 

In this work, we consider a different - and somewhat more exotic - possibility, namely, that 
the electroweak physics at a TeV is embedded in a large-extra-dimensions model of the kind 
introduced by Arkani-Hamed, Dimopoulos and Dvali (ADD) in 1998[4J. This model would 
then consist of some simple extension (with dileptonic resonances) of the Standard Model 
confined on the 'brane' - a four-dimensional hypersurface embedded in a 4+d dimensional 
'bulk' - with massless gravitons free to propagate in the bulk. On the brane, where laboratory 
experiments must be conducted, the gravitons will appear as multitudes of massive spin-2 



Kaluza-Klein states, each coupling very weakly to matter and gauge fields, but collectively 
building up to near-electroweak strengths in e + e~ collisions. We focus on the possibility of 
radiative processes, where each emitted graviton escapes detection (because of its extremely 
weak coupling to matter) and leads to missing energy and momentum signals. In this paper, 
to be specific, we have considered the possibility that an electroweak resonance, like an extra 
Z' boson, for example, could be excited in a simple process like 

e + e~ — > [x + n~ t$T 

where the missing energy (ffr) due to escaping gravitons provides the spread in the effective 
centre-of-mass energy necessary to home in on the Z' resonance. However, it may be pointed 
out at this very stage that a dileptonic resonance need not necessarily be an extra Z' boson 
- it could very well be a scalar, a tensor or an exotic spin state which are readily obtained in 
Tev-scale string theories^. 

Irrespective of the specific electroweak model with a Z' boson (provided the Z' couples to e + e~ 
and there are seven diagram topologies contributing, at the lowest order, to the process 

e + e _ — > fi + fi^G n , with exchange of, respectively, a photon, a Z boson or a Z' boson, i.e. 21 
Feynman diagrams in all. These follow from the basic s-channel boson-exchange diagram 
shown in Fig. 1, with the possibility of a graviton state being radiated, not just from an initial 
state (as shown), but from any one of the five internal or external legs or from any one of 
the two vertices. One can safely neglect the corresponding diagrams where a graviton (or a 
tower of gravitons) is exchanged in place of the neutral gauge bosons, since the corresponding 
amplitudes will be very small 1 . 




Figure 1. Feynman diagrams contributing to the process e + e — > G, 



It is necessary to make a coherent sum of all 21 diagrams to obtain the cross-section 

cr n (M n ) = cr(e + e" -» fi + fi~G n ) 

1 Snnnrr'ssorl bv four extra, nnwfirs of Mo. the strinff scale, which will he a. TeV or more. 



with radiation of a particular graviton Kaluza-Klein mode G n . The Feynman rules are given 
in Refs. jHJ U\ ■ The entire calculation, though straightforward, is long and tedious. We have 
performed these calculations by evaluating all the trace relations using the software FORMjHj. 
Each graviton contributes separately to the missing energy signal, which means that we require 
to make an incoherent sum over all graviton states which are kinematically accessible, i.e. 



making the quasi-continuum approximation for the closely-spaced graviton states. The density- 
of-states function p(M) is given in Ref. [7j. The actual cross-sections are calculated by in- 
serting the squared, spin-summed matrix elements into a Monte Carlo event generator, which 
provides numerical results of the required degree of accuracy. In this article we present only 
the physics results. 

In our numerical analysis, we have implemented the following kinematic cuts: 

1. The final state muons should not be too close to the beam pipe, which is ensured by 
demanding that the muon pseudorapidity should satisfy 77^ < 2.5 for both the muon 



2. The muons should also have significant transverse momentum, which is ensured by 
demanding that pj, > 10 GeV for both muon tracks. 

3. There should be substantial missing energy, which is ensured by demanding that $ T > 
20 GeV. It is also useful to impose an upper cut tf T < 200 GeV to get rid of photon- 
exchange processes where the muons are very soft. 

4. The dimuon invariant mass = +p A1 -) 2 , which should show a peak at the Z' 
resonance, should satisfy M MA , > 200 GeV which not only removes the peaks due to 
the photon and Z-exchange effects, but also reduces the SM background 2 substantially. 
It may be noted that we would not miss out on resonances in the 100-200 GeV range 
because these are already ruled out by LEP-2 data. 

For the graviton masses and couplings, we have adopted the minimal ADD model, which is 
characterised by two parameters, viz. the number of compact dimensions d and the (so-called) 
string scale Ms- The constraints on Ms from purely terrestrial experiments [§] vary between 
about 0.7 — 1.2 TeV, depending on d. We have, in general, chosen Ms > 1 TeV. In this paper, 
2 Here 'SM background' refers to the SM embedded in an ADD-type scenario. 
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in keeping with the general philosophy of laboratory-based analyses, we choose to ignore the 
strong astrophysical bounds[10] on the d = 2 scenario. 



There is a much wider choice when it comes to deciding the parameters of the Z' sector. In 
general, we can write the couplings of the Z' boson to fermions (/) as 



A?/ 7 = */7" [gT/PL + gT/P R \ (i) 

where Pl/r = §(1 =F 7s) and g'^ and g^}' '. The couplings of the Z' boson to chiral fermions 
are different for each scenario in) considered. In order to have a focussed discussion, we have 
chosen to follow the pragmatic approach of Dittmar et al ^1] which is limited to just five 
scenarios, of which the first three are unrealistic and the last two are well motivated. 

1. The couplings of the Z' boson are identical with those of the Z boson; the only difference 
lies in the mass M' z . The £ + £~Z' vertex (where £ = e,/i) is given by g'j^ ~ —0.20 and 

g'S - o.i7. 

2. The coupling of the Z' boson is purely vectorlike, with the vector coupling equal to 
that of the Z boson; the mass M' z is again a free parameter. The £ + £~ Z' vertex has 
gfi ^ -0.20 and g'^ = 0. 

3. The coupling of the Z' boson is purely axial, with the axial coupling equal to that of the 
Z boson; as before, the mass M' z is a free parameter. The £ + £~ Z' vertex has g'j^ = 
and g'fM ~ 0.17. 

4. The Z 1 boson arises in a SU(2)l x U(1)y x U(1)y> model which originates in the spon- 
taneous breakdown of an E 6 group: 

Eq 
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i 
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SU{3) C x SU{2) L x U(1) Y x U(1) Y > 
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SU(3) C x SU(2) L x U(1) Y 
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SU(3) C x U(l) 

em 



In this model, there are two Z'-bosons. We focus on the lighter Z\ which is a mix- 
ture Z' = Z' x cosf3 + Z'^smft, where (3 varies between — n/2 to +ir/2, with (3 = 
— tan -1 ^5/3 ~ —0.91 corresponding to direct breaking of Eq in superstring-inspired 
models, with no intermediate SO (10) group. The £ + £~Z' vertex now has g'^f ~ 
0.238 cos(P - 0.406) and g'$ ~ -0.119 sin(/3 - 0.659). 

5. The Z' boson arises in a left-right symmetric model based on a gauged SU(2) L x 
SU(2) R x U(1)b-l symmetry, with a much simpler breakdown chain 

SU(3) C x SU(2) L x SU(2) R x U(1) B -l 

i 

SU(3) C x SU(2) L x U(l) Y 

i 

SU(3) C x 17(1) 

em 

which is presumably the end sector of the breakdown of a bigger symmetry which unifies 
colour interactions with the electroweak sector. In this model, the £ + £~ Z' vertex has 
g'jj) ~ and g'jy) ~ (1 — a 2 LR ). where the parameter cxlr lies between ^2/3 ~ 
0.8165 and ^2 ~ 1.4142. 

In terms of these parameters, we can write the partial decay width of the Z' boson to a pair 
of fermions as 



( 9 2?) 2 + W) 2 -x {( 9 2;>) 2 + ( g £>) 2 - e« 



0(l-4r) (2) 



where x = (m.f/M' z ) 2 and iV c is the number of colours of the fermion /. Since, in general, 
M' z > 100 GeV from the LEP-2 constraints, we can neglect x for all but the top quark and 
get a simplified form 

\ I'., r / ., . , ■> 

(3) 



ff c 24tt 



The detailed couplings of the Z' to fermions in these five scenarios are given in Table 1. 
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Table 1. Couplings of the Z' boson to different fermions in the five models under considera- 
tion. We write cw = cos 6 W and sw = sin 9 W . 



It is now vital to numerically evaluate the total decay width of the Z' boson in these models 
and determine if, indeed, it will appear as a narrow resonance. This is illustrated in Figure 2, 
where we have plotted the total decay width Tz> as a function of M' z for different values of 
the couplings in the two realistic models in question, i.e. (4) and (5). The hatched area 
corresponds to the E 6 model, i.e. (4), while the dotted area corresponds to the left-right 
symmetric model, i.e. (5). These areas have been generated by varying the couplings f3 and 
ctLR over the full allowed range (see above). The three solid lines show the width in the 
(unrealistic) comparison models (1), (2) and (3). Small deviations from linearity in all the 
curves correspond to opening-up of the tt decay channel. The dark-shaded area is ruled out 
by the direct LEP-2 constraints. 
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Figure 2. Decay width of the Z' boson in five different models as a function of its mass. Numbers in 
parentheses correspond to the number of the model as described in the text. The dark-shaded region is ruled 
out by LEP-2. 

It may be seen that in the two realistic models (4) and (5), the decay width of the Z' boson does 
not rise above 17 GeV, even when the mass reaches 1 TeV. Since this indicates Tz>/M z > < 2%, 
we may safely conclude that the Z'-boson will appear as a narrow resonance in the discussion 
that follows. 

In any of the above models, the probability of exciting a Z' resonance in the process e + e~ — > 
fi + fi~G n will depend on the probability of emitting initial-state gravitons carrying away the 
requisite energy from the electron-positron pair. This means that the size of the resonant 
cross-section will depend rather sensitively on the ADD parameters, viz., d, the number of 
compact dimensions, and M S) the string scale. This is illustrated in Figure 3, where we plot 
the total cross-section as a function of M S) for a fixed set of all the other parameters, mostly 
chosen to get a high cross-section. Models 2 and 3 yield cross-sections in the same ballpark 
as Model 1, and have been omitted to avoid making the figure clumsy. 
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Figure 3. Cross-section for e + e~ — > o& a function of Ms, keeping d = 2 and Mz> — 250 GeV. 

The black curves correspond to Model 1, i.e. the couplings of the Z' boson to fermions are identical to the 
couplings of the Z boson. Blue curves correspond to the E^-based Model 4 with [3 — 0, the last being chosen 
to correspond to large cross-section. Red curves correspond to a LR- symmetric Model 5 with cx^n — 
which is also chosen to get the largest cross- section. 

It immediately becomes clear from Figure 3 that the present process is viable at a 500 GeV 
collider, with an integrated luminosity of around 1000 fb , only if we assume that d = 2 and 
Ms is not much greater than 1.8 GeV. For d — 3, one would require a much higher luminosity 
of at least 10 6 fb _1 in order to have observable effects. Higher values of d are simply not 
accessible at the planned energies and luminosities. 

The above curves correspond to a Z'-boson mass of 250 GeV. For the same choice of mass, 
and with d = 2 and Ms = 1 TeV, we have plotted, in Figure 4, the distribution in dimuon 
invariant mass M w = — p^-) 2 for all the five models in question. 
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Figure 4. Bin-wise distribution of the differential cross-section in dimuon invariant mass in different 
models. The dashed histogram corresponds, in every case, to absence of any Z' resonance, i.e. to a SM on 
the brane. The top row illustrates the cases for Models 1, 2 and 3. The middle row illustrates Model 4 for 
different choices of mixing angle (3, while the bottom row similarly illustrates the case of different choices of 

&LR- 



In Figure 4, the dashed histogram represents the case when there is no Z' resonance in the 
given range, i.e. one observes only the SM on the brane. This falls off sharply as M MAt increases 
because of the usual s-channel suppression far away from the photon and Z resonances. The 
presence of a Z' boson changes this behaviour dramatically, as shown by the solid histograms. 



The first three, corresponding to Models 1, 2 and 3 respectively, have large and wide reso- 
nances, consistent with the larger width and couplings illustrated in Table 1 and Figure 2. The 
behaviour of 'SM-like' and 'axial' couplings for the Z' boson is very similar simply because 
the vector coupling of charged leptons to the Z boson is very small in the SM. There is a 
qualitative difference between these and a purely vectorial coupling, in that the latter does not 
show any destructive interference effects. More realistic results are shown in the blue and red 
histograms. The blue histograms correspond to the Eq model with choices of (3 = 0, 7r/6 and 
— 7r/6 respectively. Though the model allows (5 to vary from — n/2 to +pi/2, we found that 
the qualitative behaviour of the invariant mass distribution shows very little variation beyond 
P = ±7r/6. It is quite clear that (3 = provides the largest resonances. For finite values of j3 
the resonances look small in the figure, but are not actually so, since the ordinate is plotted 
on a logarithmic scale. Similar arguments hold for the red histograms in the bottom row, 
which correspond to Model 5, with different values of the coupling constant a LR = 0.85, 1.15 
and 1.40 respectively, spanning the allowed range. As expected from the 1 /aLR dependence 
shown in Table 1, the largest resonances correspond to the smallest value of aLR- 

We reiterate at this point that the actual point in the parameter space of the model incorpo- 
rating Z' bosons is a matter of detail and our work does not focus on this. Figure 4 simply 
illustrates the fact that with quite realistic choices of parameters in this sector, we can expect 
large resonances in the dimuon mass distribution in the process e + e~ — > J^t, provided 

the model is embedded in an ADD-like brane world scenario. 

In all of the above analyses, the mass of the Z' boson was set to 250 GeV. It is natural to 
ask whether it would be possible to distinguish resonances when M z > is larger and approaches 
the machine energy of 500 GeV. Obviously, in this case one would get smaller and broader 
resonances, which may not show such spectacular deviations from the SM expectations as 
shown in Figure 4. In this case, it is better to have a quantitative index, for which we have 
adopted a x 2 fit to the line-shape. This has been done as follows: the differential cross- 
section has been calculated in 14 bins of 20 GeV each, starting from M w = 200 GeV to 
= 480 GeV. Denoting the cross-section in the ith bin by Oi and choosing a luminosity 
L = 10 3 fb" 1 , we then calculate the number of events expected in each bin as iVj = La-i. 
We separately calculate the expected number for N- SM ^ and Nf', i.e. without and with the 
Z' resonance. If the expected number drops below unity, we set it to unity, to take care of 
random fluctuations and the fact that the law of large numbers is clearly not valid. The x 2 
is now calculated, assuming Gaussian errors (statistical only) 6N$ S M) = \J N\ SM \ using the 



simple formula 

14 (n z '-N (sm) ) 2 

The value of this \ 2 is a measure of the deviation of the line-shape in M w from the SM 
expectations and consequently falls as the resonances become smaller. For x 2 < 23.7, the 
line-shape is consistent with Gaussian random fluctuations in the SM at 95% confidence level 
(C.L.), which means that this value of x 2 corresponds to the 95% discovery limit of the 
resonance. 
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Figure 5. 95% C.L. discovery limits in the Mg-Mz' plane corresponding to three different models. The 
black curve corresponds to a Z 1 boson with SM-like couplings to fermions. The blue curve corresponds to an 
Eq -based model with j3 = and the red curve corresponds to a LR symmetric model with cxlr — y/2/3. 

In Figure 5, we have plotted these 95% C.L. discovery limits in the M$-Mz> plane, the region 
below the curve being accessible to a 500 GeV linear collider with 10 3 fb _1 luminosity. The 
three curves all correspond to d — 2 and, respectively, Model 1 (black), Model 4 (blue) with 
(3 — and Model 5 (red) with a^ = 0.85, these being (as we have seen) the choices which 
lead to the largest cross-sections. As expected from the large cross-sections in Figure 3, the 
discovery limits are the most robust in the (unphysical) case of Model 1, and the ones for the 
more realistic cases look disappointing by contrast. Nevertheless, it should be possible, for 
a string scale not much larger than a TeV, to detect Z 1 resonances all the way up to nearly 
330 GeV. Given the kinematic constraints at a machine energy of 500 GeV, this is quite a 
reasonable effect. Moreover, as stated above, the nnrnose of this work is to hiffhlieht the 




possibility of exciting resonances of perhaps unknown nature, and not to focus on the details 
of models with Z' bosons. 

In this paper, therefore, we have investigated the possibility that, at a 500 GeV e + e~ collider, 
invisible graviton radiation in a brane world scenario could provide a useful tool to discover 
low-lying resonances coupling to e + e~ and /U + //~ pairs. Not only are such resonances predicted 
in extensions of the standard electroweak model, but some might be low-lying stringy excita- 
tions arising in a theory with TeV strings. We have explored the collider phenomenology of 
this process by focussing on extra Z' bosons and showed that reasonably optimistic discovery 
limits are predicted in two of the more popular models. In so doing, we have also established 
the methodology for investigation of this particular signal. It may be noted that at a 500 GeV 
e + e~ collider, one of the easiest processes to look for will be a hard muon pair with substantial 
missing energy, and we shall surely have data on this when the machine actually runs. The 
possibility that a new resonant state will be discovered — or, at least, confirmed — in this 
process is by no means a far-fetched one, and this is the substance of the present work. Quite 
obviously, there are severe limitations due to the limited energy and luminosity available at 
the machine in question, and one could carry out a more extensive study if, for example, 
the collisions took place at y/s = 2 TeV and an integrated luminosity of 10 6 fb -1 could be 
achieved |12j. However, we feel that it is premature to carry out any numerical analysis based 
on such machine parameters. Our point is sufficiently conveyed by the rather spectacular res- 
onant structures shown in, for example, Figure 4, and a study of the remaining aspects could 
await the development of a detailed machine design for the International Linear Collider. 
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